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摘要 :以 黄土 沟 密 区 典型 小 流域 泥 河 沟 为 研究 区 ,基于 1986 年 的 彩 红外 航空 相片 .2002 年 SPOT 影像 .2016 年 GF-1 卫星 影像 的 
解 译 结 果 和 社会 经 济 统计 数据 ,利用 景观 指数 土地 利用 程度 ,信息 人 等 方法 分 析 了 该 流域 近 '30 年 的 景观 格局 演变 规律 ,并 借 
鉴 生态 服务 价值 当量 估算 法 定量 探讨 了 生态 系统 服务 价值 的 变化 特征 。 结 果 表 明 : 中 近 30 年 来 研究 区 土地 利用 景观 格局 发 生 
了 显著 变化 ,除了 林地 和 建设 用 地 景观 面积 增加 ,耕地 、 园 地、 未 利用 地 均 有 不 同比 例 的 减少 名流 域内 景观 整体 破碎 度 减少 ， 
优势 斑 块 的 连通 性 呈 增 加 趋势 。@) 土 地 利用 程度 整体 呈 上 升 趋势 有 高 于 全 国平 均 水 平 231, 证 地 利用 信息 粹 先 减少 后 增加 , 流 
域 景观 经 历 了 “无 序 - 有 序 -无 序 ” 的 变化 。@ 研 究 期 内 流域 的 总 生态 系统 服务 价值 呈 持 续 上 升 趋势 ,单项 生态 服务 功能 主要 为 
土壤 形成 与 保护 .废物 处 理 .水源 涵 养 和 生物 多 样 性 保护 。 高 分 辨 率 卫 星 影 像 为 流域 景观 格局 演变 和 生态 服务 价值 分 析 提 供 了 
较为 详细 的 数据 支撑 。 
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Abstract: A case study was carried ‘out-in Nihegou, a typical watershed in the Loess Gully region. Color-infrared aerial 
photographs taken in 1986; SPOT images created in 2002, results interpreting GF-1( Gaofen- 1 satellite ) images in 2016, 
and socio-economic data were collected and used in this study. The change in landscape patterns over the past three decades 
was characterized by ‘the landscape index and land use level indices and analyzed using information entropy theory. The 
characteristics of the~change in ecosystem service values are quantitatively assessed with a version of an equivalence 
estimation method, which was originally used for ecology service valuations. The results showed that over the past 30 years, 
land use and land cover in the studied areas have undergone significant changes, with the areas of arable land, orchards, 
and tnused land decreasing and the area of woodland and developed land increasing. The overall fragmentation of the 
landscape pattérns is decreasing and the connectivity of the dominant patches is tending to grow stronger. The level of land 
useris continuing to grow and is greater than the national average level, at 231. The information entropy of land use 
decreased and then increased, indicating that the landscape patterns changed from an originally disordered state to a more 
ordered state and then returned to a more disordered state. The aggregate value of ecosystem services in this area, with the 


major individual functions being soil formation and protection, waste disposal, conservation of water sources, and 
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maintenance of bio-diversity, is tending to rise. High-resolution satellite images provide highly specific data in support of 


this characterization of the changes in the landscape patterns and the ecosystem service values of this watershed. 


Key Words: land use; landscape pattern; landscape index; ecological service value; equivalence method; Nihegou 


river basin 


景观 格局 演变 主要 关注 景观 的 结构 和 功能 随时 间 变 化 的 过 程 与 规律 ,研究 景观 格局 演变 是 探讨 生态 状 
况 、 空 间 变异 规律 以 及 与 生态 过 程 相关 的 区 域 资 源 环境 问题 的 重要 途径 。 生 态 系统 服务 功能 指 生态 系统 
与 生态 过 程 形成 及 维持 人 类 生存 发 展 的 自然 环境 条 件 '” ,对 其 价值 进行 定量 评估 是 高 效 、 合 理 配置 资源 环境 
的 基础 ” ,也 是 制定 生态 环境 建设 相关 政策 的 必要 条 件 ”。 社 会 经 济 发 展 人 类 活动 和 政策 等 因素 改变 了 景 
观 格局 中 的 物质 循环 能量 流动 和 土壤 水 分 等 '”" 生态 过 程 ,并 对 生态 系统 服务 价值 (ESV) 产 生 了 重要 影响 。 
目前 ,对 景观 格局 的 研究 大 多 集中 在 景观 格局 变化 及 驱动 力 分 析 方 面 "” ,对 生态 系统 服务 的 研究 主要 是 关 
注 ESV 的 空间 变化 “" ,或 者 是 侧重 于 湿地 农田 森林“" 等 某 一 生态 系统 服 ,将 景观 格局 演变 和 ESV 响 
应 进行 结合 的 研究 还 有 待 深入 。 

黄土 高 原 是 地 球 上 分 布 最 集中 且 面 积 最 大 的 黄土 区 ,同时 也 是 世界 上 水 土 流失 最 严重 和 生态 环境 最 脆弱 
的 地 区 之 一 。 为 了 防治 水 土 流失 ,恢复 生态 平衡 建立 良性 循环 农业 生态 系统 ;国家 在 “七 五 "时 期 将 “黄土 高 
原 综合 治理 研究 ”纳入 重点 科研 项 目 “ 区 域 综合 治理 试验 ”之 中 ;并 在 黄 士 高 原 不 同类 型 区 设置 了 11 个 试验 
示范 区 。 泥 河沟 流域 作为 黄土 沟 包 区 的 典型 代表 ,经 过 “七 五 "一 “ 九 五 ”的 连续 玫 关 ,其 景观 格局 和 农户 的 经 
济 生活 状况 发 生 了 巨大 变化 ,有 必要 进行 长 期 的 追踪 探讨 ,但 目前 对 该 流域 的 研究 多 集中 在 水 土 流失 治理 、 粮 
食 作物 增产 .土壤 水 肥 条 件 改善 "| 等 方面 ,尤其 是 “ 九 五 "重点 攻关 结束 以 来 ,由 于 国家 对 该 流域 的 治理 力 
度 有 所 下 降 , 导 致 与 之 相关 的 研究 也 相应 减少 ,本文 以 泥 河沟 流域 为 研究 对 象 ,将 景观 格局 演变 和 ESV 响应 
相 结 合 来 揭示 黄土 沟 密 区 这 种 特殊 地 貌 近 30 年 来 ESV 的 演变 规律 ,这 对 深入 认识 科技 攻关 成 果 及 后 续 资源 
的 合理 配置 ,实现 黄土 高 原 经 济 与 生态 环境 的 协调 发 展 具 有 重要 意义 ,同时 为 国家 今后 在 黄土 高 原 地 区 挖掘 
资源 潜力 .改善 生态 条 件 .繁荣 区 域 经 济 和 解决 水 土 流 失 等 方面 制定 政策 提供 理论 参考 。 


1 研究 区 概况 


泥 河沟 流域 (图 1) 位 于 黄土 高 原 南部 陕西 省 咸阳 市 淳 化 县 境内 , 属 典 型 的 黄土 沟 蜜 区 。 流 域 面积 9.48 
km ,包括 七 个 行政 村 ,2016 年 常住 人 口 约 为 3607 人 。 地 势 北 高 南 低 ,海拔 699 一 1163m , 沟 殉 上 游 坡 陡 沟 罕 ， 
侵蚀 活跃 ,下游 侵蚀 趋 缓 ; 年 平均 气温 9.8% ,无 霜 期 183d, 太阳 年 辐射 总 量 504.35KJ/cm , 为 暖 温 带 半 湿 润 
气候 。 年 均 雨 量 600:6mm,7 一 9 月 份 降 水 量 占 全 年 的 50% 以 上 且 多 为 暴雨 。 自 然 植 被 多 为 草本 灌木 日 位 于 
沟 坡 。 主 要 景观 类 型 为 耕地 .园地 林地 .草地 水域 .建设 用 地 和 未 利用 地 。 


2 数据 来 源 与 方法 


2.1 数据 源 及 预 处 理 

遥感 影像 分 别 是 1986 年 彩 红 外 航空 相片 (1:10000) .2002 年 2.5m 空间 分 辨 率 的 SPOT 影像 和 2016 年 的 
GF-1 卫星 影像 (2m 分 辨 率 全 色 /8m 分 辨 率 多 光谱 ) ,利用 ENVI 5.3 对 SPOT-5 和 GF-1 数据 进行 辐射 定 标 、 
FLAASH 大 气 校正 、 正 射 校正 .融合 以 及 和 1986 年 由 彩 红外 航 片 解 译 的 土地 利用 现状 图 进行 配 准 , 配 准 误差 
控制 在 0.5 个 像 元 以 内 ,然后 利用 面向 对 象 分 类 方法 对 两 景 高 分 辩 率 影像 进行 人 机 交互 解 译 ,根据 全 国 农 业 
区 划 委 员 会 1984 年 发 布 的 《土地 利用 现状 调查 技术 规程 》 制 定 的 土地 利用 现状 分 类 及 含义 并 结合 流域 特点 ， 
将 土地 利用 类 型 划分 为 7 种 ,分 别 为 耕地 ( 旱 耕 地 和 轮 荡 地) .园地 (果园 ) .林地 (乔木 林 、 疏 林地 和 造林 地 ) 、 
草地 ( 天然 草地 和 人 工 草地 ) 水域 (水 库 .渠道 等 ) .建设 用 地 (农村 居民 点 、 独 立 企业 用 地 和 道路 ) 和 未 利用 
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地 (陡坡 荒草 地 ,上 废 沟 滩 地 等 ) 。 基 于 混淆 矩阵 方法 和 
实地 验证 进行 解 译 精度 评价 ,3 期 影像 解 译 精度 均 大 于 
90% ,满足 分 析 的 需要 。 此 外 还 包括 淳 化 统计 年 鉴 和 实 
地 走访 数据 等 数据 。 
2.2 ”研究 方法 
2.2.1 景观 格局 演变 

(1) 景 观 特征 指数 选取 ”景观 指数 指 高 度 浓缩 景 
观 格局 信息 ,反映 其 结构 组 成 和 空间 配置 某 些 方面 特征 
的 简单 定量 指标 。 根 据 需要 在 景观 类 型 水 平 上 选取 
斑 块 类 型 面积 (CA) 、 斑 块 面积 比例 (PLAND) 、 斑 块 数 
(NP) 景观 形状 指数 (LST) .聚集 度 (CLUMPY) 、 斑 块 密 
度 (PD) .最 大 斑 块 指数 (LPI) .平均 斑 块 面积 [AREA_ ER 
MN) ;景观 水 平 上 选取 斑 块 数 (NP) 、 斑 块 密度 (PD) 香 ,| 水系 
农 多 样 性 (SHDI) .香农 均匀 度 (SHEI) 景观 形状 指数 a 
(LSI) 蔓延 度 (CONTAG) .景观 丰富 度 (PR) .平均 斑 块 i 6, 
面积 ( AREA_MN) .最 大 斑 块 指数 (LPI) ,以 上 指标 可 以 0 04 08miles 
反映 各 景观 类 型 面积 .破碎 度 .形状 复杂 度 .均衡 性 等 ， 
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二 
计算 及 意义 详 见 r91 。 3 
这 oO 了 9 o: 了 O r 
(2) 土地 利用 动态 度 ”单一 土地 利用 动态 度 指 某 “六 0 ee 本 
研究 区 一 定时 间 范围 内 某 种 土地 利用 类 型 的 年 变 a 
化 率 '”" o Fig.1 Sketch map of the location of the study area 
L, a m 1 
De x— x 100% (1) 
L 了 


式 中 :D 为 研究 时 段 内 某 一 土地 利用 类 型 的 单一 动态 度 ;L, .L, 分 别 为 起 始 年 和 终止 年 某 景观 类 型 的 面积 ;7 
为 研究 时 上段。 
(3) 土 地 利用 程度 ”土地 利用 程度 用 于 分 析 区 域 土 地 利用 的 广度 和 深度 ,能 够 定量 化 反映 自身 自然 属性 
和 人 类 与 自然 环境 的 综合 效应 "1 
L, = 100 x YP x 0, (2) 
式 中 .LL 为 土地 利用 程度 综合 指数 ;P, 为 第 i 级 土地 利用 程度 分 级 指数 ;0, 为 第 i 级 土地 利用 程度 分 级 面积 
分 比 。 
(4) 二 地 利用 信息 烂 “土地 利用 信息 入 可 综合 反映 某 区 域 在 一 定时 段 内 各 土地 利用 类 型 的 动态 变化 及 
其 转换 程度 ,对 于 区 域 土地 利用 结构 调整 具有 指导 意义 , 炉 值 越 大 ,土地 利用 系统 的 有 序 程度 越 低 ,各 要 素 发 
展 越 均衡 ,反之 亦 然 , 当 系统 完全 达到 平衡 时 ,信息 炉 达到 极 大 值 '*”，。 
H=- > Pilogp, (3) 
式 中 ;P, 为 研究 区 各 景观 面积 占 该 流域 总 面积 的 百分比 ;n 为 景观 类 型 数量 。 
2.2.2 ESYV 评估 
目前 ,国内 外 对 ESV 评估 的 方法 主要 以 Costanza 等 5 中 提出 的 生态 价值 估算 方法 为 主 , 谢 高 地 等 25 在 
Costanza 评估 模型 基础 上 ,制定 了 中 国生 态 系 统 生态 服务 价值 当量 因子 表 '* 。ESYV 当量 因子 的 价值 量 估算 公 
式 为 : 
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1 二 
rc = 开 xPx 一 2 0， (4) 
了 Nn i=1 


式 中 :VC 为 ESV 当量 因子 的 价值 量 (元 km “a ') ;P 为 平均 粮食 价格 (元 /kg) ;0 为 研究 区 平均 粮食 产量 
(kg/km ) ;n 为 年 份 数 。 本 文 把 耕地 和 园地 统一 对 应 于 生态 资产 类 型 中 的 农田 ,以 研究 区 平均 粮食 产量 
281930kg/km? 作为 基准 单产 ,并 结合 平均 粮食 价格 1.68 元 /kg, 考 虑 到 在 没有 人 力 投 入 的 自然 生态 系统 提供 
的 经 济 价值 和 现 有 单位 面积 农田 提供 的 食物 生产 服务 价值 是 1:7 的 关系 ,得 出 泥 河 沟 农田 自然 粮食 产量 的 经 
济 价值 VC 为 67663.2 元 kg ” a  。 价 值 当量 乘 以 流域 农田 自然 粮食 产量 的 经 济 价值 即 可 定量 化 出 泥 河沟 景 
观 类 型 单位 面积 生态 价值 。 评 估 泥 河沟 的 ESV 的 公式 为 : 


ESV = > (A. x VC.,) (5) 


ESV, = 2 (4; x VCn) (6) 
式 中 :ESYV 是 研究 区 生态 系统 服务 价值 (元 ) ;4, 为 景观 类 型 上 的 面积 (km?) ,VC 为 -ESV 系数 (元 km ”a '); 
ESV ,为 生态 系统 第 了 项 服务 功能 价值 ;VC 是 景观 类 型 的 第 f 项 服务 功能 价值 系数 (元 km ”a ')。 
表 1 泥 河沟 流域 各 地 类 生态 系统 单位 面积 服务 价值 系数 /( 元 km-? a-0) 


Table 1 Ecological service value coefficients per unit area of different landscape types ‘in Nihegou 


景观 类 型 Tandscapefpe 


生态 系统 服务 
Reosystem services 耕地 园地 林地 草地 水 域 未 利用 地 
Arable land Orchard Woodland Grassland Water body Unused land 

气体 交换 Gas exchange 33831.60 33831.60 236821.20 54130.56 0.00 0.00 
气候 调节 Climate regulation 60220.25 60220.25 182690.64 60896.88 31125.07 0.00 
水 源 涵养 Water conservation 40597.92 40597:92 216522.24 54130.56 1378976.02 2029.90 
0 os 本 98788.27 98788.27 263886.48 131943.24 676.63 1353.26 
废物 处 理 Waste treatment 110967.65 110967-65 88638.79 88638.79 1230116.98 676.63 
生物 多 样 性 保护 Biodiversity conservation 48040.87 48040.87 220582.03 73752.89 168481.37 23005.49 
食物 生产 Food production 67663.20 67663.20 6766.32 20298.96 6766.32 676.63 
原材料 Raw material 6766;32 6766.32 175924.32 3383.16 676.63 0.00 
娱乐 休闲 Entertainment and leisure 676.63 676.63 86608.90 2706.53 293658.29 676.63 
合计 Total 467552.71 467552.71 1478440.92 489881.57 3110477.30 28418.54 


2.2.3 ”敏感 性 验证 


为 验证 选用 的 TC 是否 适合 本 研究 区 ,选用 敏感 性 指数 ( CS) 来 分 析 。 敏 感性 分 析 指 在 求 得 某 个 模型 的 最 
优 解 后 ,研究 模型 中 某 个 或 若干 个 参数 允许 变化 到 多 大 , 仍 能 使 原 最 优 解 的 条 件 保 持 不 变 !”] 。 通 过 计算 CS 
来 确定 EESV 随时 间 变 化 对 VC 的 依赖 程度 ,CS 指 VC 变动 1% 引 起 ESV 的 变化 状况 ,CS>1, 表 明 ESV 对 VC 是 
富有 弹性 的 ,结果 的 准确 性 越 低 ;相反 则 表明 ESV 对 VC 是 缺乏 弹性 的 55 ,结果 可 信 。 此 外 ,即使 CS<1, VC 
的 值 过 高 或 过 低 都 会 一 定 程度 上 影响 ESV 随时 间 变 化 的 准确 性 。 现 将 生态 服务 价值 系数 分 别 调整 50% 来 计 
算 -05. 


_ | (ESV, = ESV,) /ESV, 

(VC = VG) AVC; 
式 中 ,C5 为 敏感 度 , 指 VC 变动 1% 引 起 的 ESV 的 变化 情况 ;VC 表示 生态 价值 系数 ;i 和 j 表示 初始 价值 系数 和 
调整 后 的 价值 系数 ;k 指 各 景观 类 型 。 


CS (7 
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3 结果 与 分 析 


3.1 景观 类 型 水 平 变化 特征 分 析 

由 表 2 看 出 泥 河 沟 流域 耕地 的 形状 指数 (LSI) 呈 先 减少 后 增加 的 趋势 , 斑 块 密度 持续 增加 ;园地 斑 块 类 型 
面积 至 2002 年 增加 幅度 达到 24.04% ,2002 一 2016 年 则 减少 27.73% , 斑 块 密度 上 升 ,聚集 度 先 是 明显 上 升 随 
后 又 呈 下 降 趋 势 , 景 观 形状 指数 总 体 呈 减少 趋势 ,表明 景观 破碎 度 减 少 ,人 类 的 开垦 利用 使 园地 的 形状 趋 朱 规 
则 , 斑 块 连通 性 增加 ;林地 的 斑 块 面积 呈 持 续 增 加 趋势 ,其 聚集 度 在 各 地 类 中 表现 最 高 ;草地 的 斑 块 数 和 形状 
指数 呈 持 续 减 少 趋势 ; 随 着 居民 生活 水 平 提高 ,房屋 .道路 和 厂房 等 建设 用 地 面积 逐步 增加 ,同时 未 利用 地 得 
到 了 有 效 的 利用 ,水 域 在 研究 期 内 变化 不 大 。 
3.2 景观 级 别 水 平 变化 特征 分 析 

从 ( 表 3) 看 出 该 流域 景观 形状 指数 (LSI) 在 1986 一 2002 年 间 呈 显著 下 降 趋势 ,说明 人 类 活动 干预 将 景观 
形状 改造 趋 于 集中 。 草 延 度 (CONTAG ) 在 研究 期 内 旦 上升 趋 势 ,景观 整体 破碎 化 程度 略 有 降低 。 香 农 多 样 性 
此 数 ( SHDI) 和 香农 均匀 度 指数 (SHET) 均 表 现 为 先 减少 再 略微 增加 的 趋势 ,香农 多 样 性 指数 (SHEI) 减少 , 表 
明 各 景观 类 型 在 景观 分 布 中 均衡 化 程度 较 低 , 而 香农 均匀 度 指数 的 减少 则 说 明 流 域内 景观 优势 度 在 增加 , 景 
观 优势 类 别 对 景观 整体 的 控制 作用 增强 ,主要 是 林地 斑 块 面积 增加 ,对 景观 的 控制 作用 增强 导致 的 。 


表 2 1986 一 2016 年 泥 河沟 流域 景观 类 型 水 平 指数 
Table 2 Indexes on landscape class level of the Nihego River watershed from 1986 to 2016 


类 型 年份。 机关 关 。 班 类 面 。 部 g、 时 观 政 于 全 度 天 夫 区 。 最 大 殖 。 区 面积 
Type Yo 型 积 积 比例 证 状 指数 人 MP Pp 块 指数 AREA_MN/ 
CA/km2 PLAND/% LSI LPL/% i 
耕地 Arable land 1986 348.7619 37.2284 28 11:4387 0.892 2.9888 16.1087 12.4558 
2002 304.0734 C32.6563 18 1.9331 13.5879 9.622 16.893 0.9109 
2016 296.4341 31:6497 42 4.4843 8.8725 10.6189 7.058 0.9008 
园地 Orchard 1986 158.3433 16:9023 35 8.933 0.9076 3.7361 5.0901 4.5241 
2002 196.4211 “21.0949 17 1.8257 5.738 6.812 11.5542 0.9361 
2016 141.9609 15.1569 39 4.164 4.2666 9.0606 3.64 0.9028 
林地 Woodland 1986 149.849 15.9955 80 15.2353 0.8308 8.5396 2.0165 1.8731 
2002 231.3709 一 一 妇 .8483 18 1.9331 18.3267 11.185 12.8539 0.8917 
2016 263.0196 28.0821 26 2.776 23.1724 10.5222 10.1161 0.9006 
草地 Grassland 1986 115.9152 12.3733 90 17.2682 0.7895 9.607 4.4406 1.2879 
2002 74.6097 8.0128 26 2.7923 3.0925 9.5347 2.8696 0.8688 
2016 103.5676 11.0577 25 2.6692 7.0655 9.9408 4.1427 0.8795 
水 域 Water body 1986 8.8428 0.9439 6 5:20 0.8176 0.6405 0.6758 1.4738 
2002 7.4118 0.796 1 0.1074 0.796 3.8261 7.4118 0.8663 
2016 8.5786 0.9159 1 0.1068 0.9159 3.3673 8.5786 0.8967 
建设 用 地 1986 48.4831 5.1753 32 8.2069 0.8656 3.4158 1.1687 1.5151 
Build-up land 2002 54.6714 5.8715 54 5.7994 1.7472 9.3065 1.0124 0.853 
2016 74.348 7.938 31 3.3098 1.4119 7.9722 2.3983 0.892 
未 利用 地 1986 106.6222 11.3813 31 16.4244 0.7937 3.3091 2.7651 3.4394 
Unused land 2002 62.5745 6.7203 41 4.4032 1.9939 11.3258 1.5262 0.8287 
2016 48.6993 5.1995 49 5.2316 1.6877 9.7241 0.9939 0.8381 


CA: Class area; PLAND : percent of landscape; NP: Number of patches; LSI: Landscape shape index; CLUMPY : Clumpiness; PD: Patch density; 
LPI: Largest patch index; AREA: Area 
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表 3 1986 一 2016 年 泥 河沟 流域 景观 级 别 水 平 指标 
Table 3 Indexes on landscape level of the Nihegou River watershed from 1986 to 2016 


年 份 下 类 禾 。 瑞 度 香农 习性 香农 均 认 pi J 2 
ear NP PD SHDI SHEI TSI CONTAG PR ARFA_MN LpI/% 
1986 302 32.2368 1.6647 0.8555 16.4126 44.4768 7.0000 3.1020 16.1087 
2002 175 18.7943 1.6284 0.8368 12.2874 48.331 7.0000 5.3208 18.3267 
2016 213 22.7416 1.6481 0.8469 12.6169 47.2997 7.0000 4.3972 23.1724 


NP: Number of patches; SHDI: Shannon's diversity index; SHEI: Shannon's evenness index; CONTAG : Contagion index; PR: Patch richness ; 
AREA_MN : Mean area 


3.3 ”景观 类 型 面积 变化 特征 及 相互 转化 关系 分 析 

图 2 为 研究 区 3 期 土地 利用 类 型 图 ,由 ( 表 4) 看 出 研究 区 各 景观 类 型 面积 均 发 生 了 变化 ,耕地 呈 持 续 减 
少 趋 势 , 变 化 量 为 -0.5377km ;园地 先 增加 后 减少 ,变化 量 为 -0.1622km ;林地 变化 较 显 著 , 从 1.5069km 增加 
到 2.6323km? ,变化 量 为 1.1254km? ;建设 用 地 占 比 从 5.14% 增 加 到 7.84% 而 未 利用 地 则 从 占 比 12.07% 减 少 到 
6.31% ;草地 和 水 域 面积 变化 不 大 。 整 体 看 , 近 30 年 间 耕 地 、 园 地 和 林地 占 优势 地 位 ;其 面积 之 和 在 3 个 时 期 
都 超过 了 研究 区 总 面积 的 70%。 由 图 3 来 看 ,林地 、 草 地 和 未 利用 地 表现 出 较 明显 的 变化 速率 ,其 中 林地 达到 
3.34%。 男 外 ,除了 林地 ,建设 用 地 动态 变化 为 正 外 (>0) ,其 他 地 类 均 为 负 动 态 变化 。 由 表 5 看 出 变化 最 明显 
的 是 林地 与 园地 和 耕地 之 间 的 相互 转化 , 占 整 个 研究 区 变化 面积 的 18.80%’ 该 转换 关系 一 方面 是 受 政策 影 


0 0.4 0.8m 
t | | 


图 例 


目 园地 国 国 水域 攻 国 林地 半天 草地 
证 未 利用 地 故国 建设 用 地 恒 辐 耕地 


图 2 泥 河 沟 流域 1986、2002 和 2016 年 土地 利用 类 型 图 
Fig.2 Land use types in Nihegou River watershed in 1986, 2002 and 2016 
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响 ,部 分 坡 耕 地 和 一 些 坡度 较 大 容易 导致 水 土 流失 的 园地 被 有 计划 地 转换 为 林地 ,恢复 林 草 植被 ,使 得 林地 
面积 增加 ; 男 一 方面 ,原来 一 些 品种 落后 的 果树 和 老 树 逐 渐 被 经 济 效益 较 好 的 新 品种 果树 和 耕地 取代 ,导致 园 
地 的 数量 和 结构 发 生变 化 ,此 外 , 随 着 泥 河 沟 治理 的 深入 推进 ,加 上 农业 技术 承包 、 良 田 改 造 等 的 实施 和 市 场 
经 济 的 影响 ,使 得 流域 内 的 景观 布局 尤其 是 耕地 园地 和 林地 的 比例 得 到 调整 和 优化 。 


表 4 各 类 景观 面积 变化 量 及 百分比 
Table 4 The landscape area change and percentage of the Nihegou River watershed 


只 佬 忆 

od 1986 年 2002 年 2016 年 
Landscape type 面积 可 放 队 2 职 可 

ea a a ee Re 人 
耕地 Arable land 3.5008 36.93 3.0382 32.05 2.9631 31.26 -0.4626 “一 00751 0.5377 
园地 Orchard 1.5829 16.70 1.9633 20.71 1.4207 14.99 0.3804 -0.5426. 门 “0.1622 
林地 Woodland 1.5069 15.90 2.3129 24.40 2.6323 27.77 0.8060 0.3194 1.1254 
草地 Grassland 1.1692 12.33 0.7477 7.89 1.0370 10:94 £0.4215 0.2893 -0.1322 
水 域 Water body 0.0890 0.94 0.0743 0.78 0.0851 090 -0.0147 0.0108 -0.0039 
建设 用 地 Build-up land ”0.4871 5.14 0.5470 5.77 0.7432 7.84 0.0599 0.1962 0.2561 
未 利用 地 Unused land 1.1441 12.07 0.7963 8.40 0.5986 631 -0:3478 -0.1977 -0.5455 
合计 Total 9.4800 1.00 9.4800 1.00 9.4800 1.00 0 0 0 


表 5 1986 一 2016 年 景观 类 型 面积 转移 矩阵 /km 
Table S Transition matrixes the Nihegou River watershed landscape area from 1986 to 2016 
2016 年 景观 类 型 Landscape type in 2016 


和 社 地 园地 林地 草地 水 域 。 ”建设 用 地 未 利用 地 加 各 
Arable land Orchard Woodland ™~ Giassland Water body Build-up land Unused land 
耕地 Arable land 2.106 0.508 0:491 0.124 0.272 0.003 0.032 1.430 
园地 Orchard 0.440 0.778 0.209 0.078 0.129 一 0.007 0.863 
林地 Woodland 0.156 0.057 0:878 0.287 0.045 0.009 0.064 0.618 
草地 Grassland 0.120 0.038 0.495 0.328 0.009 0.005 0.162 0.829 
水 域 Water body 0.104 01034 0.045 0.022 0.281 一 0.000 0.205 
建设 用 地 Build-up land 0.004 0.000 0.024 0.002 一 0.045 0.013 0.043 
未 利用 地 Unused land 0.032 0.003 0.557 0.251 0.005 0.021 0.207 0.869 
转 入 面积 Transfer-in area 0.856 0.640 1.821 0.764 0.460 0.038 0.278 4.857 


3.4 土地 利用 程度 分 析 

泥 河 沟 3 期 土地 利用 程度 总 体 呈 上 升 趋势 且 高 于 全 国平 均 水 平 231( 表 6)。 说 明 “ 七 五 ”一 “ 九 五 ”时 期 
的 连续 攻关 ;使 该 流域 的 土地 利用 方式 趋 于 合理 ,土地 利用 程度 大 大 提高 。 后 15 年 呈 略 微 下 降 趋 势 是 由 于 随 
着 经 济 发 展 和 城市 化 进程 加 快 ,外 出 务工 人 员 增 加 ,人 们 对 土地 缺乏 管理 ,使 土地 没有 得 到 有 效 利用 ,此 外 , 建 
设 用 地 一 定 程度 上 破坏 了 景观 的 完整 性 和 连通 性 ,景观 功能 有 所 下 降 。 


表 6 泥 河 沟 流域 1986 一 2016 年 土地 利用 程度 综合 指数 
Table6 The integrative indexes of land use degree from 1986 to 2016 in Nihegou River watershed 
年 份 Year 1986 2002 2016 
土地 利用 程度 指数 Land use degree 251.86 255.90 255.64 


3.5 土地 利用 信息 坑 评价 
由 公式 (3) 得 出 3 期 土地 利用 信息 箭 分 别 为 0.7248 .0.7148 .0.7219 ,特征 为 先 减少 后 增加 且 总 体 略 有 下 
降 ,1986 年 的 土地 利用 信息 炉 最 高 ,表明 土地 利用 系统 有 序 度 较 低 ;2002 年 最 低 则 表明 土地 利用 系统 的 有 序 
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度 增加 ,无 序 度 降低 。 信 息 灶 的 变化 表明 泥 河沟 土地 利 0 

用 系统 经 历 了 “无 序 一 有 序 一 无 序 ” 的 逐步 调整 过 程 。 | 

经 过 “七 五 "时 期 的 集中 治理 ,该 流域 的 土地 利用 结构 总 | 

得 到 调整 ,使 得 2002 年 的 土地 利用 系统 有 序 性 与 治理 。 襄 ， 

前 相 比 明显 提高 ,土地 利用 系统 逐渐 稳定 。 二 江 加 

3.6 生态 系统 服务 价值 演变 分 析 3 | 置 民 SN N 

3.6.1 总 生态 系统 服务 价值 变化 由 a 
近 30 年 来 , 泥 河 沟 流域 总 ESV 呈 持 续 上 升 趋势 ， 全 “| 二 

由 548.69x104 元 增加 到 673.11x104 元 ,增加 了 22.68% 0 

( 表 7)。 从 各 景观 类 型 的 ESV 来 看 ,1986 一 2002 年 间 ， 她 一 民 

除 林地 和 园地 的 ESV 呈 增 加 趋势 外 ,其 他 景观 类 型 的 i 

ESV 均 呈 减少 趋势 ,其 中 林地 的 ESV 增加 119.16x10” 图 3 1986 一 2016 年 土地 利用 动态 度 


元 增 幅 达 到 53.49% ;2002 一 2016 年 间 林地 草地 和 水 Fig.3 Dynamic change of land use from 1986 to 2016 
域 的 ESV 均 呈 增加 趋势 。 总 体 上 只 有 林地 的 ESV 保持 
增加 趋势 , 且 近 30 年 总 增长 率 达 到 74.68% ,林地 单项 ESV 的 增加 对 该 流域 总 ESV 的 增加 起 重要 作用 。 


表 7 泥 河 沟 流域 不 同类 型 景观 的 生态 系统 服务 价值 


Table 7 Different landscape types ecosystem services value of Nihegou River watershed 


4 TT/s 年 各 人 
a ESV/x104/( 元 /a) 1986 一 2002 2002 一 2016 年 1986 一 2016 年 
恋 仁 最 恋 仙 最 恋 仙 县 
Landscape type 1986 年 。 2002 年 。 2016 年 。 肥 化 于 “ 变 优 率 / 史 一 受 化 一 变化 率 /% 和 化 最 训 化 率 /% 
(x104 元 ) (x104 元 ) (x104 元 ) 
耕地 Arable land 163.68 142.05 138.54 二 21.63 -1321 -3.51 -2.47 -25.14 -15.36 
园地 Orchard 74.01 91.79 66.43 17. 访 24.02 -25.37 -27.64 -7.58 —10.24 
林地 Woodland 222.79 341.95 389.17 119.16 53.49 #7 3 13.81 166.38 74.68 
草地 Grassland 57.28 36.63 50.80 =20:65 -36.05 14.17 38.68 一 6.48 -11.31 
水 域 Water body 27.68 23.11 26.47 -4.57 -16.51 3.36 14.54 -1.21 -4.37 
利用 
未 利用 地 3.25 2.26 1:70 -0.99 -30.46 -0.56 -24.78 -1.55 -47.69 
Unused land 
总 计 Total 548.69 637:79 673.11 89.11 16.24 35.31 5.54 124.42 22.68 


3.6.2 单项 生态 系统 服务 价值 变化 

泥 河沟 流域 单项 ESV 变化 较 明显 ( 表 8) ,除了 食物 生产 的 价值 呈 下 降 趋势 ,其 他 单项 ESV 均 有 不 同比 例 
的 增加 。1986 一 2002 年 间 , 原 材料 和 娱乐 休闲 的 ESV 变化 较 其 他 功能 大 ,变化 率 分 别 为 46.06% 、39.09%、 
27.92% ;与 前 15 ,年 相 比 ,2002 一 2016 年 间 , 食 物 生 产 、 气体 交换 水源 涵 养 ESV 变化 较 显 著 。 从 泥 河 沟 流域 
ESV 的 构成 来 看 ,各 单项 功能 价值 所 占 比例 的 大 小 依次 为 :土壤 形成 与 保护 > 废物 处 理 > 水 源 涵养 > 生物 多 样 
性 保护 > 气候 调节 > 气体 交换 > 原材料 > 食物 生产 > 娱乐 休闲 。 泥 河沟 土壤 形成 与 保护 的 功能 价值 最 为 突出 , 占 
总 功能 价值 的 19.86% ;废物 处 理 . 水 源 涵养 及 生物 多 样 性 保护 的 功能 价值 次 之 ,三 者 各 占 总 功能 价值 的 13% 
左右 ;娱乐 休闲 占 总 功能 价值 的 比例 最 低 ,为 3% 左右。2015 年 位 于 泥 河 沟 咀 头 村 的 现代 农业 示范 园区 “首届 
葡萄 文化 节 ” 开 幕 ,该 园区 以 休闲 农业 旅游 为 载体 , 以 观赏 品尝 .采摘 日 光 温室 大 棚 新 品种 葡萄 ( 京 蜜 、 京 香 
玉 ` 红 巴 拉 多 等 ) 为 重点 ,吸引 了 大 量 游 客 来 体验 现代 新 型 农业 科技 管理 等 文化 ,促进 了 该 流域 娱乐 休闲 功能 
的 逐步 提高 。 
3.6.3 ”敏感 度 分 析 

由 (图 4) 看 出 ,CS 从 低 到 高 依次 是 未 利用 地 ,水域 .草地 .园地 耕地、 林地 且 均 小 于 1。 其 中 ,林地 的 敏感 
性 指数 最 高 ,2016 年 为 0.58 ,说 明 当 林地 的 生态 价值 系数 增加 1% 时 ,生态 系统 服务 总 价值 将 增加 0.58% 。 敏 


http ://www.ecologica.cn 


ChinaXiv 合 作 期 刊 


8 期 宋 敏 敏 ” 等 :黄土 沟 窄 区 小 流域 景观 格局 演变 及 生态 服务 价值 响应 2657 


感性 分 析 表 明 , 泥 河沟 ESV 对 于 价值 系数 变化 缺乏 弹性 ,生态 价值 系数 可 靠 ,研究 结果 可 信 ,可 为 该 流域 的 生 
态 建设 提供 参考 。 


表 8 泥 河 沟 流域 1986 一 2016 年 ESV 结构 变化 


Table 8 Changes of structure of ecosystem services value of Nihegou River watershed during 1986 一 2010 


ESVZ(x104 元 ) 1986 一 2002 2002 一 2016 1986 一 2016 
景观 类 型 
变化 量 / 变化 量 / 变化 量 / 
Landscape type 2 恋 化 玄 / 护 恋 化 谈 / 号 恋 化 谈 7] 纺 
1986 2002 2016 (xl04 元 ， 变化 率 /% (xl104 元 ) 变化 率 /% (x104 元 ) 变化 率 /% 
与 休 交 
飞 体 交换 59.21 75.74 82.78 16.53 27.92 7.04 9.29 23:37 39.81 
Gas exchange 
与 候 调 节 
飞 候 调节 65.54 77.16 81.07 11.62 17,73 3.91 5.07 15.53 23. 
Climate regulation 
水 源 涵养 3 72.1 84.84 92.26 12.74 17.67 7.42 8.75 20.16 27.96 
Water conservation 
土壤 形成 与 保护 
Soil formation and 105.57 120.42 126.54 14.85 14.07 6:12 5.08 20.97 19.86 
protection 
废物 处 理 91.16 91.82 91.68 0.66 0.72 —0.14 -0.15 0.52 0.57 
Waste treatment 
生物 多 样 性 保护 70.42 83.64 89.58 13.22 18.77 5.94 71 19.16 27.21 
Biodiversity conservation 
食物 生产 ， 37.93 37.03 33.65 -0.9 F237 -BB8 -9.13 一 4.28 —11.28 
Food production 
原材料 . 30.35 44.33 49.63 13:98 46.06 5.3 11.96 19.28 63.53 
Raw material 
娱乐 休闲 16.4 22.81 25.91 6.41 39.:09 3.1 13.59 9.51 57.99 
Entertainment and leisure 
合计 Total 548.68 637.79 673.1 89.11 16.24 35.31 5.54 124.42 22.68 
疆 愉 已 二: 0.7 
4 结论 与 讨论 
四 1986 
> 口 || 旦 、 和 已 0 国 2002 
ESV 估算 大 致 分 两 类 ,分 别 是 基于 单位 服务 功能 价 9 0 2016 
y 机 = 号 0.5 
格 法 和 基于 单位 面积 价值 当量 因子 法 '*” ,不 同方 法 
选取 的 模型 和 参数 都 会 使 结果 存在 差异 ,目前 还 没有 让 加 04 
公众 普遍 接受 的 ESV 评估 体系 。 考 虑 到 前 者 输入 参数 缠 芋 03 
较 多 .计算 过 程 复杂 且 对 每 种 服务 价值 的 参数 标准 难以 2， 
统计 ,而 后 者 较为 直观 县 适用 于 区 域 ESV 的 评估 i ， 加 
因此 本 文选 用 后 者 ,在 谢 高 地 等 (= 加 的 生态 服务 价值 :三国 
1 广 当 yA ge RE 0 人 Ee Zh 
当量 估算 法 的 基础 上 ,探讨 1986 一 2016 年 泥 河 沟 流 域 耕地 ”园地 林地 草地 ”水 域 未 利用 地 
ESV 的 变化 特征 ;同时 利用 景观 指数 、 土 地 利用 动态 度 景观 类 型 Landscape type 
中 言 息 ， 4 和 泪 定 量 4 其 景观 局 演变 见 律 得 出 以 
和 信息 依 等 方法 定量 分 析 其 景观 格局 演变 规律 ,得 出 的 站 
结论 如 下 : Fig.4 Ch f coeffici itivi i 
g. anges of coefficients of sensitivity of ecosystem services 
GD) 近 30 年 来 泥 河沟 土地 利用 景观 格局 变化 显 value from 1986 to 2016 


著 , 除 林地 和 建设 用 地 面积 增加 外 ,耕地 、 园 地、 未 利用 

地 面积 均 有 不 同比 例 的 减少 。 由 景观 指数 分 析 结 果 可 知 流域 内 景观 整体 破碎 度 减少 ,优势 斑 块 的 连通 性 呈 增 
加 趋势 。 耕 地 园地 与 林地 之 间 的 转换 是 主要 的 用 地 类 型 转化 方式 。 除 林地 、 建 设 用 地 为 正 动态 变化 外 ,其 他 
均 为 负 动 态 变 化 。 国 家 对 该 流域 的 治理 使 其 利用 结构 得 到 优化 ,土地 利用 程度 大 大 提高 且 高 于 全 国平 均 水 平 
231 ,土地 利用 系统 趋 于 稳定 有 序 。 
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(2) 从 总 ESV 变化 角度 看 , 近 30 年 的 ESV 呈 持 续 上 升 趋势 。 这 与 李 唱 和 任 志 远 2 的 陕 北 黄土 高 原 土地 
利用 ESV 时 空 研究 和 李娜 等 :5 探讨 的 黄土 丘陵 区 ESV 的 结果 一 致 。 从 各 景观 类 型 的 ESV 来 看 ,1986 一 2002 
年 间 ,林地 和 园地 的 ESV 呈 增 加 趋势 ,其 他 景观 类 型 的 ESV 均 呈 减少 趋势 ,2002 一 2016 年 间 , 林 地 .草地 和 水 
域 的 ESV 均 呈 增加 趋势 。 从 单项 ESV 变化 角度 看 ,1986 一 2016 年 间 该 流域 单项 ESV 变化 较 明 显 ,除了 食物 
生产 的 价值 呈 下 降 趋势 外 ,其 他 单项 ESV 均 有 不 同比 例 的 增加 。 从 流域 ESV 的 构成 看 ,各 单项 功能 价值 中 以 
土壤 形成 与 保护 .废物 处 理 水源 涵 养 和 生物 多 样 性 保护 为 主 , 其 中 土壤 形成 与 保护 的 功能 价值 最 为 突出 ; 占 
总 功能 价值 的 19.86%。 经 分 析 , 泥 河沟 ESV 变化 的 原因 主要 是 景观 格局 导致 各 自然 生态 系统 面积 发 生变 化 
和 国家 政策 的 推动 。 

(3) 通 过 敏感 性 指数 验证 ,该 流域 ESV 对 生态 服务 功能 指数 缺乏 弹性 ,生态 价值 系数 的 不 确定 性 对 整个 
研究 区 ESV 估算 结果 的 稳定 性 影响 不 大 ,研究 结果 可 靠 ,能 够 为 深入 认识 科技 攻关 成 果 及 后 续 资 源 配 置 提供 
参考 。 然 而 ,要 想 准 确 获 取 该 区 域 ESV ,唯一 方法 就 是 采用 特定 区 域 生 态 系统 服务 的 特定 单位 :22 ,因此 本 
研究 仍 需 深入 研究 。 
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